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Q ■ Decoherence in Nature has become one of the most pressing problems in 

physics. Many appHcations, including quantum information processing, depend 
on understanding it; and fundamental theories going beyond quantum mechan- 

l|-l3| , where the breakdown of quantum theory appears 



^ . ics have been suggested 

^ \ as an 'intrinsic decoherence', mimicking environmental decoherence [4]. Such 

t^^ ■ theories cannot be tested until we have a handle on ordinary environmental 

o: 

decoherence processes. Here we show that the theory for insulating electronic 
spin systems can make accurate predictions for environmental decoherence in 
^ ■ molecular-based quantum magnets [5|]. Experimental understanding of decoher- 

H : . 

c^ ' ence in molecular magnets has been limited by short decoherence times, which 

nn 

make coherent spin manipulation extremely difficult [6H9|]. Here we reduce the 
decoherence by applying a strong magnetic field. The theory predicts the con- 
tributions to the decoherence from phonons, nuclear spins, and intermolecular 
dipolar interactions, for a single crystal of the Fcg molecular magnet. In ex- 
periments we find that the decoherence time varies strongly as a function of 
temperature and magnetic field. The theoretical predictions are fully verified 



experimentally - there are no other visible decoherence sources. Our investiga- 
tion suggests that the decoherence time is ultimately limited by nuclear spins, 
and can be extended up to about 500 microseconds, by optimizing the temper- 
ature, magnetic field, and nuclear isotopic concentrations. 

Environmental decoherence processes are reasonably well understood at the atomic scale 
10| (although some poorly understood noisy sources remain [11|). However both quantum 



information processing, and the fundamental tests noted above, require an understanding 
of decoherence in larger systems, where experimental decoherence rates are usually much 
larger than theoretical predictions. This discrepancy is usually attributed to "extrinsic" 
sources (external noise, uncontrolled disorder /imp urities) . We thus need to find systems, 
with many degrees of freedom, where extrinsic decoherence can be eliminated, and where 
we have a quantitative understanding of other decoherence sources. 

Many insulating electronic spin systems are currently the subject of intense experimental 
interest, notably in semiconductor quantum dots [I4I, nitrogen- vacancy centres in diamond 
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-Il5|. and large spin magnetic molecules [6 



8[. In all these systems, three environmental 



decoherence mechanisms are involved. The electronic spins couple locally to (i) phonons (an 



oscillator bath 16|); (ii) to large numbers of nuclear spins (a spin bath (l7|); and (iii) to each 



other via dipolar interactions. The long range of dipolar interactions is a major problem: it 



makes quantum error correction more difficult, is theoretically complicated 18|], and is very 
hard to eliminate experimentally. 



Two great advantages of the Fes system 19[ are that (i) the interaction strengths are 



well known, allowing quantitative predictions; and (ii) it can be prepared with little disorder 
and few impurities, reducing the danger of extrinsic decoherence. The number of relevant 
environmental degrees of freedom is very large; depending on isotopic concentrations, there 
are 10^° — 10^^ nuclear spin levels in each molecule; and the system couples to a bulk phonon 
bath. At low temperature, the electronic spin system behaves as a 2-level quantum bit 
(qubit), with a splitting depending strongly on an applied transverse field, perpendicular to 
the easy axis z (See Fig. 1(b)). While almost all previous experiments in the electron spin 
systems examined the low-field regime, where nuclear spins make the qubit dynamics inco- 
herent, here we go to high transverse fields; by suppressing most nuclear spin decoherence, 
we unveil the coherent qubit dynamics. 



Spin echo measurements, using a Hahn echo sequence 20[, were carried out on two single 



crystals of Feg molecules (see Fig. 1), with natural isotopic concentrations, using a 240 GHz 
pulsed ESR spectrometer [2l|, |22]. Fig. 1 shows Feg ESR signals obtained by echo-detected 
field sweep measurements. The resonant peaks are broadened, with a width ~ 0.1 T; the 
peculiar structure of the peak when H || y comes from dipolar interactions 23|. 
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FIG. 1: Typical spin-echo detected field-sweep ESR signals as a function of transverse magnetic 
field H_|_ recorded in two different samples, at different temperatures and orientations of field, (a) 
For sample 1, Solid red line - H_|_||y, T = 1.58 K. Dashed blue line - H_|_||i;, T = 1.67 K. The sample 
dimensions are approximately z : x : y = 1000 : 700 : 250 ^iv? . The lower left inset (i) shows the 
low-T spin structure of the Feg molecule. The lower right inset (ii) shows the directions of the easy 
and hard axes of a Feg molecule (a approximately gives the direction of the crystallographic vector 
a), (b) For sample 2, Solid red line - H_|_||y, T = 1.7 K. Dashed blue line - H_|_||x, T = 1.23 K. 
The sample dimensions are approximately z : x : y = 900 : 800 : 400 iim^ . The lower inset shows 
the tunneling splitting, 2Ao, as a function of transverse field at H_|_||y (solid red line) and at H_|_||x 
(dashed blue line). 



These ESR signals may be understood as follows. The qubit splitting 2Ao varies with field 



as shown in the inset to Fig. 1(b). For fields Hy = 9.5T, iJ^. = 0, or H^ = 11.3T,Hy 



0, 



the electronic Hamiltonian for Fes 



23| predicts 2Ao(Hj 



11.5 K, equivalent to our 



spectrometer frequency of 240 GHz, implying we should see resonance peaks at these fields. 
These predictions are reasonably well satisfied in both samples. The discrepancies, discussed 



in detail in the supplementary material, come from two sources: (i) sample misorientation 



(ii) weak high field departures from the model Hamiltonian 2J] used to predict the field 
splitting. 

In this experiment, the uniform ESR precession mode (a k = magnon) interacts with its 
surroundings; to understand the decoherence time T2, and the ESR lineshapes, we need to 
look at the processes contributing to them. For convenience, we define a dimensionless deco- 
herence rate 7^ = hj I\^i- Then the contributions to 7,^ are as follows (the full quantitative 
discussion, for the 2 samples in this experiment, is in the supplementary material): 

(i) In this high field regime, nuclear spins interact locally with each molecular spin, and 
cause decoherence by a 'motional narrowing' process in which they attempt to entangle with 



the fast-moving qubit [18|, l25[. The nuclear decoherence rate is 7^^ = £'q/2A^, where Eo is 
the half- width of the Gaussian multiplet of nuclear spin states coupled to the qubit; and the 
nuclear contribution to the ESR linewidth is just E^. Now in this experiment, with naturally 
occurring isotopic concentrations, i?o = 4 x 10~^ K^ whereas Aq = 5.75 K. Thus 7^ ~ 10~^ 
is quite negligible, and the nuclear spin contribution ~ Eq to the linewidth is also very small 
compared to the main contributions. Isotopic substitution of deuterium for the 120 protons 
in each molecule will further decrease 7^ by a factor 15.2 to 7^ ~ 7 x lO"'^^. 

(ii) The form of the local spin-phonon interaction is determined by the system symmetry. 
At high fields this interaction simplifies, and we find a dimensionless phonon decoherence 
rate 18[ given by 7^ = [(J-A5A^)/(7rpc^^'^)] coth(Ao/A;T), where p is the sample density, c^ 



the sound velocity, and Tas is the spin part of the spin-phonon matrix element 23| . While 
7? depends strongly on field, its temperature dependence is very weak for these experiments. 
The contribution of this spin-phonon process to the ESR linewidth is negligible. 

(iii) The intermolecular dipole interaction leads to direct coupling of the k = ESR 
precession mode to finite momentum magnons; these cause strong decoherence. The fc = 
magnon can in principle decay spontaneously into multiple magnons, or scatter off ther- 
mal magnons. The dipolar interaction affects the ESR lineshape and the decoherence rate 
very differently. The long-range interaction creates a distribution of demagnetization fields 
around the sample. In the case of highly polarized samples, this is strongly sample-shape 
dependent, but for annealed samples, it is Gaussian distributed J26|, |27|; in both cases it 
can be calculated numerically. The lineshape then refiects the quite broad distribution of 
these fields. However, decoherence comes from the magnon decay process described above. 



and depends only on the phase space available at the resonance field for these processes; 
it can then be calculated directly from the analytic expression for this process. At the ex- 
perimental temperatures, the magnon decoherence rate ~ exp{—2Ao/kBT), coming almost 
entirely from thermal magnon scattering. This decoherence contribution increases rapidly 
with temperature. 

(iv) Finally, there can be extrinsic contributions from impurities and defects, which 
typically cause the easy-axis anisotropy parameter D of the Feg Hamiltonian to fiuctuate 



around the sample (the '/^-strain' effect 28|, |29|). This will certainly contribute to the ESR 



linewidth. However these impurities and defects will not contribute to the decoherence rate 
at all, provided they are static, since then they simply shift the individual qubit energies. 
Any impurities or defects with significant dynamics will cause extrinsic decoherence, which 
would then lead to a disagreement between theory and experiment - which is not seen in 
the present study. 

The resulting theoretical decoherence times for Fes, for each sample, are presented in 
Fig. 2, together with the results of the measurements. The agreement is very good; we 
emphasize that apart from the size of the spin-phonon coupling, which is only known ap- 
proximately, there are no adjustable parameters in these fits. 

A number of features should be stressed here. First, notice how different the effects of 
these mechanisms is on decoherence and on the the ESR lineshape. The ESR linewidth and 
lineshape are completely dominated by static impurity fields and by the spatially varying 
dipolar fields. However the decoherence is completely dominated, at these high fields, by the 
phonons and dipolar interaction effects. At lower fields, the nuclear spin decoherence would 
also be important, but its effect on the ESR linewidth would still be negligible. Note also 
that whereas the dipolar contributions to the ESR lineshape depends strongly on sample 
shape, this shape can only affect the dipolar decoherence near the edges of this line. In 
the middle of the line, where the decoherence is coming from molecules near the centre 
of the sample, where the field is homogeneous, one expects no dependence of the dipolar 
decoherence on sample shape. This is also what we found in the experiment. 

Second, we emphasize how the experiment tests the phonon and dipolar contributions 
to the decoherence separately - they have very different temperature dependences in the 
regime covered here, with phonons dominating below ~ 1.2 K, and magnons dominating at 
higher temperature. We find agreement between theory and experiment, with no adjustable 
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FIG. 2: The measured and calculated decoherence times T2 in samples 1 and 2, as a function of 
temperature, (a) Results for H||y. Thin line with diamonds - experiment in sample 1, Hy = 9.845 
T. Thin line with circles - experiment in sample 2, Hy = 9.875 T. Vertical and horizontal bars 
represent the experimental errors and uncertainty in temperature (AT = ±0.05 K). Thick solid 
line - theory including phonon and magnon contributions, Hy = 9.5 T. The inset shows the partial 
contributions calculated for T2 (solid line) from magnons (dashed line) and phonons (long-dashed 
line), together with the corresponding experimental results for the two samples (diamonds and 
circles). The scale on the right-hand side of the main figure indicates the decoherence Q-factor 
Qfj, = vr/70 = 7rT2Ao/ti; the right-side scale on the inset shows ^^. (b) The same as for (a), but now 
for H||x. The experimental curves were measured at H^ = 10.865 T (sample 1) and H^ = 11.953 
T (sample 2). The theoretical curves are obtained at H^ = 11.3 T. Again, the scales on the right- 
hand side of the figure and inset indicate the decoherence Q-factor Q^p and 7^ respectively, as in 
(a). 

parameters, across nearly an order of magnitude in decoherence rate. Thus all decoherence 
in the experiment can be accounted for by environmental sources. This implies "we have no 
measurable extrinsic decoherence here, either from disorder or dynamic impurities. Nor do 
we have evidence for any other contributions, either from "3rd-party decoherence" 4], or 
from any of the "intrinsic decoherence" sources [l| discussed in the literature. 

The present investigation suggests that T2 in the Fcg SMM system can be increased 
significantly by optimizing the temperature and the qubit splitting 2Ao. If we extend the 
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FIG. 3: Calculated contributions to the decoherence coming from the coupling to nuclear spins, 
phonons, and magnons respectively, (a) shows the 3 individual contributions which sum to give the 
dimensionless decoherence rate 7,^ = h/T2Ao, as a function of the qubit splitting in the case H_|_| |rE. 
(b) shows the corresponding three contributions to the decoherence time T2. In each figure: Brown 
dashed line - nuclear contribution; Solid lines - magnon contributions at different temperatures 
from 0.1 — 1.6 K; Long dashed lines - phonon contributions, shown for the same temperatures as 
in the magnon case. 



calculations to cover all temperature and field regimes accessible to resonance experiments, 
we obtain 23| the results shown in Fig. 3. Going to lower temperature allows use of a 



smaller ESR frequency 2Ao; both these changes strongly reduce the dipolar and phonon 
decoherence contributions, giving a large increase in T2. The optimal decoherence rate 
comes at the point where the phonon and nuclear spin decoherence contributions cross (for 
2Ao ~ 2 K); provided that kT < 0.13 K, so that the dipolar contribution can be ignored. 
One then finds that 7^ ~ 1.5 x 10^^, corresponding to a decoherence time T2 ~ SO/xsecs. 
Any further increase in T2 would require isotopic substitution, notably of deuterons in place 
of the protons; we then find the optimal decoherence time to be T2 ~ 500/isecs, at 2Ao = 
0.8K = 17 GHz, and T = 45 mK. These considerations show the usefulness of this kind of 
theory in the optimal design of spin qubit systems. If quantum mechanics is to be tested 
on anything but microscopic scales, it will be essential to continue developing theory and 



8 

experiment for systems like the present one, where the environmental decoherence processes 
can be understood quantitatively, and where extrinsic decoherence sources can be largely 
eliminated. 
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SUPPORTING ONLINE MATERIAL 

Decoherence in Quantum Magnets: Methods 

S. Takahashi, I. S. Tupitsyn, J. van Tol, C C Beedle, D. N. Hendrickson and P. C E. Stamp 

o' 

Here we provide details of the methods used in our high-frequency pulsed ESR experi- 
ments, and in the theoretical analysis of environmental decoherence contributions. 

in 

^ : (i) EXPERIMENTAL METHODS: We investigated two slab-shaped single crystals 

of Fcg magnetic molecules (SMMs). These Peg crystals were synthesized using the original 

^ . method developed by Weighardt et al. [l|. ESR measurements were carried out using a 240 

GHz continuous- wave (cw) /pulsed ESR setup built at the National High Magnetic Pield 
Laboratory (NHMPL), in Tallahassee PL, USA J2|, yj. The system consists of a 40 mW cw 
solid-state source at 240 GHz, quasioptical ESR bridge, superheterodyne detection system 
and 12.5 tesla superconducting magnet, and employs a "^He flow cryostat to vary sample 
temperature down to ~ 1.2 K. Por control of the sample orientation relative to the external 
magnetic fields, the sample crystal was placed in a rotating sample holder consisting of a 
teflon rotating rod and nonmagnetic gears. The sample angle can be changed by rotating a 

t^^ , control rod connecting from room temperature to the gear in the sample holder |2| . 

O I To measure the spin decoherence time T2, we employed a Hahn echo sequence (7r/2 — r — 

t^^ ■ TT — r — echo) where the delay time r is varied j^. 7r/2 and n are square pulses whose width 

(typically 200 ns to 300 ns long) was optimized for the maximum echo signals. This width 
corresponds to approximately 0.5 Gauss of the excitation bandwidth, much smaller than the 



c^ 



o 



o 



^ . ESR linewidth of the lowest lying Peg transition. The decoherence time T2 was extracted 

c^ \ from the decay rate of the echo area, which was well fit by a single exponential exp(— 2r/T2). 

The temperature dependence of T2 at frequency of 240 GHz was then measured between 1.2— 
2K, with the external magnetic field aligned either along the hard (H||x) or intermediate 
(H||y) crystal axes. Above 2 K, T2 became too short to give spin echoes within the time 
resolution of the pulsed ESR spectrometer. 

Por the angle- dependent ESR measurement, the sample crystal was carefully mounted 
in a rotating sample holder, allowing us to vary the direction of the external magnetic field 



along the xy-plane. The orientation of the crystallographic axes was initially identified 
through X-ray diffraction measurements, and then the low-T determination of the actual 
sample orientation was made using cw ESR measurements at 20 K. To accurately determine 
the sample orientation, we rotated the sample crystal through 360° relative to the direction 
of the magnetic fields, and fit the measured resonance fields to the calculated values as a 
function of angle with the Feg spin Hamiltonian given in eqtn. ([1]) below. We show here 
an example of such a fit, for sample 2 (see Fig. 1). As a result of these fits, we determined 
that the direction of the external magnetic field in the sample 1 was tilted by < 14° from 
the xy-plenae, whereas that in sample 2 was tilted by only < 3°. 
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FIG. 1: Plot of the resonance fields measured for Sample 2 as a function of rotation angle of the 
sample holder. The dashed red line shows the values calculated from the Feg spin Hamiltonian in 
eqtn. ([T]), described herein. 



(ii) THEORETICAL METHODS: The Feg molecule has eight Fe^+ ions Sj, of spin 
5/2 each, with an electronic spin Hilbert space of dimension 6^ ~ 1.7 x 10^. However below 
about 10 K these are locked by superexchange inte^ractions into a 'giant spin' array, of total 
spin |S| = 10. This has the effective Hamiltonian 



nis) = -DSl + ESl + Ki{SX + St) - 7S ■ H_ 



(1) 



where the easy/hard axes are along z and x respectively. The Zeeman coupling 7 = gefJ'BfJ'O, 
where (?£ ~ 2 and the transverse field H^ is perpendicular to z. The parameters of ([I]) were 



tuned [6| to accurately reproduce both the experimental value of the tunneling gap 2Ao(Hx) 
and its period of oscillations with field |7|, and are D/k^ = 0.23 K, E/k-Q = 0.094 K, and 
Ki/kB = -3.28 X 10-5 K. 

The Hamiltonian ([T]) is easily diagonalized numerically. For all fields used in this exper- 
iment, a low-energy doublet clearly separates off from the higher-energy states - this is the 
electronic 'qubit'. The energy splitting of this doublet was shown in Fig. 1(b) of the main 
text. We write these 2 states as \S), \A), where the lowest energy state is \S). Let us define 
linear combinations \2±) = {\S) ± \A))/2^^'^. In the experiment, almost all molecular spins 
are in the low-T \S) state. A 7r/2 microwave pulse then sends these into the \Z+) state, 
thereby initiating spin oscillations between \Z±) at a frequency 2Ao/h (which may then 
resonate with the applied ESR signal at frequency Qr = 240 GHz, provided 2Ao/h = ^r). 
If we ignore any field inhomogeneity (which, because of dipolar interactions, does exist in 
the sample - see below), these spin oscillations are equivalent to a uniform spin precession 
around the sample along the field directions, i.e., to a /c = magnon. A calculation of en- 
vironmental decoherence is then equivalent to a calculation of the coupling of this magnon 
to the various environmental modes. 

Decoherence Rates: We now calculate the three different sources of environmental 
decoherence in the experiment (and their rather different effect on the ESR linewidth). 
Note that the results of these calculations are shown graphically in Fig. 3 of the main text. 

Nuclear Spins: In this high field regime, the dimensionless nuclear decoherence rate 
reduces to 7^ = -E'^/2A^, where -E-o is the half-width of the Gaussian multiplet of nuclear spin 
states coupled to the qubit states |6| . To evaluate Eo one requires all the hyperfine coup 
These were previously determined for the Feg system, both in low-field experiments 
and theoretically J6|. We extend the theory to high fields by defining the vector ujk, for the 
k-th nuclear spin in the molecule, with components 



mgs. 

nfi 



H 



^" = ^E^S(K)'-K)'') (2) 



2 



where {s")^"^ is the expectation value of the individual spin s^, when the molecule is in the 
\S), \A) states. Then 2a;fc measures the difference ccf — uj-^ between the hyperfine fields 
acting on this nuclear spin when the molecule is in the |iS), |^) states. The nuclear multiplet 



halfwidth is then 

k 

The hyperfine couphngs {^;^«} themselves are essentially independent of field, but the field- 
dependence of the qubit states leads to a strong decrease of Eo at high fields. This is easily 
found numerically: for these calculations we assumed natural isotopic concentrations, and 
found Eo = Ax 10~^ K for the fields used in this experiment. This gives a quite negligible 
contribution to the decoherence rate. It also gives a very small contribution Eo to the ESR 
linewidth. 

If we do selective isotopic substitution to remove nuclear spins, we can further decrease 
the nuclear contribution to the decoherence rate. Thus substitution of deuterium for protons 
decreases Eo from the above value, computed assuming protons at the 120 different Hydrogen 
sites, by a factor 3.9, further decreasing the decoherence rate by a factor 15.2. 

Phonons: The large tunneling gap means that decoherence from low-energy acoustic 
phonons is important. We define e^js = {dj^Ua + 9Q,n^)/2 and Uajs = {djjUa — daUp)/2, 
the symmetric and antisymmetric parts of the strain tensor, written in terms of the lattice 
displacement u(r) for a molecule at position r. For Feg crystals, the mass density p = 
1920kgm~^ for the natural isotopic concentrations used here, and the unit cell volume is 
Vo = 1956A . We assume a phonon spectrum Ug = qCg, with the sound velocity Cg = 
(/cB6z)/^)(K)/67r^)^/'^ assumed the same for transverse and longitudinal phonons; here 0^ = 
33 K is the Debye temperature [6|. 

We write the spin-phonon coupling as Vs,ph = J2i ViOfOf , where the subscript i labels 
different irreducible group representations for the lattice and molecular symmetry groups, 
and the operators Of , 0\ operate separately on the spin and phonon degrees of freedom 
[lO|. At high fields the dominant terms in Vs^ph can be written as 

+ {'n2.(^xy + Vi^xy){Sx, Sy} (4) 

+ {'n3(^xy + 'n^^xy){Sx, Sy} (5) 

where the indices (X), (Y) indicate that H is parallel to x or to y respectively and the 
anticommutators {A, B} = AB + BA are evaluated at the fields used here. The coupling 



constants ?2i can only be determined accurately from experiment. However, simple estima- 



tions 
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L|-|13| shows that they are ~ -D/2, where D is the main easy-axis anisotropy energy. 
Since the phonons themselves are only weakly affected by the field, all field effects come 
from the field dependence of the qubit states and AJii±). The spin-phonon decoherence 
rate is then calculated using standard methods [IJ, [l5| to give 

A3 



'ph 
J^AS = 



J=. 



AS 



npclh'^ 



coth{Ao/kBT) 



dft. 



TT 



J2v^f^m{•A\^f\S)\' 



(6) 
(7) 



where dQ is an infinitesimal solid angle in real space, and the /i(fi) are rather complicated 
angular functions coming from the tensors eap and Ua/3, which are independent of q. If we 
now go to the high field regime, and assume that all the couplings are equal to D/2, we find 



,{X) 



L>2 



n^s' « ^{\{A\{S^,S^}\S)\' + \{A\{S^,Sy}\S)\') 



riY) 



D' 



J^'as' ^ ^i\{A\{Sy,S,}\S)\'+\{A\{S^,Sy}\S)\') 



(8) 
(9) 



These functions are evaluated numerically, by exact diagonalisation; the results are shown 
here in Fig. O 



200 



§ 

c 
o 

■c 

9- 

c 

I 



100 



1 1 1 1 1 


1 1 1 1 1 


- 


Hj^=Hx 

1 j \ 




- 


1 \ \ 
1 / ^ \ 
' / ^ \ 

' / 
/ / 




■ 




■ 


1 / 
■ / / 

ii 


Hz=0 


■ 


ii 
1/ 
/ ■ , , , 1 


1 i 


- 



5 10 15 

H±(T) 

FIG. 2: Plots of Tas defined in Eqs. dH]),®, as a function of the transverse applied field. The 
results are shown for H_|_ along x (red) and along y (blue). 



In general the phonon decoherence rate depends on temperature, but in our experiments. 
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ksT ^ 2Ao, and so the rate has already reached its low-T saturation value; this is seen for 
both field orientations in Fig. 2 of the main text. 

We note immediately that the contribution hTpi^ to the ESR linewidth will be completely 
negligible. 

Dipolar Interactions: The dipole field of each molecular spin leads to a global sample 
shape- dependent demagnetization field iidm{T^), which varies around the sample. Thus the 
qubit energy splitting 2Ao(H(r)) varies around the sample with the in-plane component of 
the total local field H(r) = Ho + Hdm(i")- Only in the close vicinity of 'resonant surfaces' in 
the sample, where 2Ao(H(r)) = hQji, can molecules be fiipped by the applied ESR signal. 
The strength of the ESR signal at a given applied field Hq is then directly proportional to 
the number of molecules in this surface, i.e., to the density of states 

AA(Ho) = fd\S{hnn-2Ao{llo + lldmir)) (10) 

The lineshape will then be complicated, with multiple peaks coming from the van Hove 
singularities arising each time a new resonant surface moves in from the sample boundary, 
as one varies the external field. By further varying this field, we sweep the resonant surfaces 
through the sample. The edges of the lineshape correspond to fields where the resonant 
surfaces are small and near the sample boundaries, in regions where the demagnetisation field 
is varying rapidly with r. In the middle of the ESR line, the resonant surfaces occupy a large 
region near the centre of the sample, where the demagnetization field is relatively uniform 
in space; there are then a large number of molecules in resonance. In this latter case, the 
approximation that the tipping pulse corresponds to a uniform (in space) precession mode, 
i.e., a. k = magnon, is a very good one. We therefore compare theoretical calculations of 
the decoherence with T2 measurements taken at the centre of the line, where the signal is 
strongest anyway. 

In the present experiment, where ksT ^ 2Ao, only 4-magnon processes where the k = 
mode scatters off a thermally excited magnon, with two magnon final states, contribute to 
the decoherence; decay into three magnons is blocked by kinematic restrictions. It then 
follows that the dipole contribution to the dimensionless decoherence rate 7^ is 

27r 






^^ = ^ 1. ir^r-^t^J ^(^0 + u;q - u^' - u;q_qO (11) 



qq' 



(4) 



where ^[?T^q] is the statistical distribution function for the magnons, T L is the 4-magnon 
scattering amphtude, and the magnon dispersion Wq depends on the apphed field and on 



the sample geometry [15[. We calculated numerically the dipolar field distributions and 
the magnon dispersion and decoherence rates for samples 1 and 2, approximating them as 
parallelipipeds with the measured aspect ratios. The Fes molecules, described by the full 



-spins Sj = 5/2 model with the known positions of each Fe ion 16|, were organized in 



triclinic arrays in the numerica 
and crystallographic axes as in 



calculation of the magnon decoherence rates, with angles 

5|. We found 



1| and directions of the anisotropy axes as in 
the width of the dipolar distributions in both samples to be ~ 0.1 T which is comparable 
with the width of the experimental ESR peaks. The magnon contribution to the decoherence 
rate for this process is shown in Figs. 2 and 3 of the main text. 

Finally, as noted in the main text, the sample certainly contains some static disorder, 
coming from dislocations and defects, as well as from other lattice strains (which can be 
induced by a strong field). The disorder will also cause variation of the Fes Hamiltonian 



parameters between molecules [17|, Il8[. The spin energy levels of each molecule are thus 
in general shifted, by different amounts, by this disorder. This will certainly influence the 
sample linewidth and lineshape, and thus give an inhomogeneous broadening contribution 
to T2. However, provided this disorder can be considered to be static, it cannot give any 
contribution to the decoherence, which is completely unaffected by static level shifts. We 
also note that the experiment is designed so that any inhomogeneous broadening is canceled 
by the rephasing vr-pulse in the Hahn echo sequence. However, any dynamic local degrees 
of freedom, such as 'loose spins' {i.e., spin impurities with energy levels in the energy range 
below ~ 10 K), or mobile charge defects, can cause decoherence. No evidence is seen in the 
experiments for such an extrinsic contribution. 
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